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I. SUMMARY 


The objective of the work conducted in this program was to define both aero- 
dynamic and manufacturing design details for an advanced 4:1 pressure ratio 
single stage centrifugal compressor at a 10 Ibm/sec flow size. The approach 
selected was to perform an exact aerodynamic scale of DDA's 404- III compressor 
from its design flow of 3.655 Ibm/sec to the required 10 lbm/sec flow size. 

Design tasks performed during this program included: 

o Aerodynamic design and analysis 
o Thermal analysis of the scaled impeller 

o Structural analysis of the scaled impeller including both static stress 
and vibration analyses. 

The results of these tasks are reviewed in Sections II and III. Section IV 
presents a detailed manufacturing definition of the impeller blading and wheel 
geometry. The manufacturing definition includes a detail drawing of the im- 
peller geometry and punched card description of both "hot" and "cold" impeller 
geometry. 
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II. AERODYNAMIC DESIGN AND ANALYSIS 


The 4:1 Rc, 404- III single stage centrifugal compressor was designed in 1975 
for use in an advanced regenerative gas turbine engine for truck/bus and power 
generation applications. The impeller design combined advanced aerodynamic 
features such as high back curvature (50°) and low blade loading with geometry 
completely compatible with production casting techniques. Goal compressor 
performance was achieved on the inital build and efficiency goals were ex- 
ceeded by over 1% after one rematch. The design point total to static effi- 
ciency was 83.3% at a point with 8% minimum surge margin. This efficiency 
level is still the best total to static efficiency demonstrated in its flow 
cl ass. 

The flow size of the 404- III compressor is 3.655 lbm/sec at a mechanical speed 
of 36015 rpm. To scale the compressor to 10 lbm/sec, a 1.6529 scale factor 
must be applied to all linear dimensions. A true scale of this compressor 
would result in a diffuser exit radius of 16.306 inches. This radius exceeded 
current NASA rig constraints and, therefore, had to be reduced to 14.3 
inches. A 90° annular bena was designed to direct the flow from the radial to 
the axial direction. 

A complete meriaional elevation of the scaled compressor is shown in Fig- 
ure 1. The defined geometry is for the "hot running" conoiton with no impel- 
ler to shroud clearance adjustment. The original 404- III compressor was test- 
ed witn a smooth approach inlet bell and rotating spinner. These contours are 
shown in Figure 1 and specifically defined in Table I. 

The impeller consists of 15 full blades and 15 splitters. "Hot" flow path 
contours, tangential thickness distributions and polar coordinate definition 
of the blade meanline are presented for the hub and shroud contours in Tables 
II and III.. Table II presents full blade definition while Table III gives the 
splitter geometry. The blade surface elements are constructed linear from hub 
to shroud along the defined quasi-normals. The coordinates given in Tables II 
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Figure 1. Meridional Flowpath for Scaled Impeller. 
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FABLE I. INLET BELL AND 

SPINNER CONTOUR COORDINATES. 

INLET BELL: CIRCULAR ARC WITH CENTER 4F C -1-057 *5.974) 

coordinates: 

AXIAL 

RADIAL 

LOCATION 

LOCATION 

-2.908 

5.974 

-2.850 

5.513 

-2.750 

5.2 25 

-2.500 

4.814 

-2.250 

4.5 58 

-2.000 

4.381 

-1.750 

4 .2 57 

-1.500 

4.177 

-1.057 

4.123 

-0.050 

4.128 

0.000 

4.132 

INLET SPINNER: 

coordinates: 

AXIAL 

R AD I AL 

LOCATION 

LOCATION 

-3.3276 

0 .0000 

-3.3111 

0 .1 759 

-3.2450 

0 .3886 

-3.1623 

0.5446 

-2. 9970 

0 .7 564 

-2.8317 

0 .9091 

-2.6665 

1.0293 

-2.3359 

1 .2093 

-2. 0053 

1.3344 

-1.6747 

1.4192 

-1.3441 

1.4711 

-1.3C0C 

1.4 757 

-1.0 COO 

1 .5000 

-0.5C0C 

1.5603 

0.0000 

1 .6314 


and III are for the impeller in the "hot and running" condition. A detailed 
description of the manufacturing or "cold" geometry is presented in Section IV. 


Using data from the BU257 rig test of the 404- III compressor as a guide to im- 
peller blockage and efficiency levels, a streamline curvature intrablade aero- 
dynamic analysis was performed for the defined scaled compressor geometry. 
The aerodynamic analysis was conducted at design speed with the following in- 
put vari ables: 
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TABLE II. SCALED 404-EII IMPELLER COORDINATES - HOT BLADE 

MEAN 3.ADE DEFINITION - FULL BLADE. 
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TABLE III* SCALED 40^-ItI IMPELLER COORDIMATES - HOT BLADE 

MEAN BLADE DEFINITION - SPLITTER* 
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o Corrected flow 
o Corrected speed 
o Inlet pressure 
o Inlet temperature 


10 Ibm/sec 
21 789 rpm 
14.7 psia 
518. 7°R 


The resulting distributions of impeller relative velocity and blade loading 
distributions are shown plotted as a function of percent meridional distance 
in Figures 2 through 7. 


The vane diffuser consists of 24 modified, two-dimensional wedge vanes with 
the leading edge located at a radius ratio of 1.0778 relative to the impeller 
exit. The diffuser entrance region is shown in Figure 8 and is centered on 
the impeller exit. The diffuser has an overall area ratio of 2.754 with a 
total divergence angle of 7.791°. The vane passage cross-section is presented 
in Figure 9 with an enlarged view of the leading edge shown in Figure 10. The 
individual vanes are constructed from straight line segments between points 1 
ano 2, Figure 9, for the pressure surface ano between points 4 and 5 for a 
portion of the suction surface. The leading edge portion of the suction sur- 
face is formed by an arc as shown in Figure 10. The suction surface arc has a 
radius of curvature of 45.233 inches and is tangent to the leading edge circle 
at points 3 ano to the straight line between points 4 and 5 at point 4. The 
diffuser exit radius is 14.30 inches and dumps directly into a 90° annular 
beno. 


The annular bend is shown in Figure 11 with detailed coordinates presented in 
Table IV. Primary considerations in the design of the annular beno were: 


o Minimize static pressure gradients at the diffuser exit plane 
o Maintain maximum flowpath radius at 16.0 inches 

To avoid large static pressure gradients at the diffuser exit, the annular 
bend was designed with a generous radius of curvature to gap ratio of 2.0. 
The area distribution shown in Figure 12 was selected to reduce velocity 
levels around the bend and, thereby, reduce total pressure losses. 
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Figure 2. Impeller Hub Relative Velocity. 









Figure 3. Impeller Mean Relative Velocity. 
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Figure 4. Impeller Tip Relative Velocity. 
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Figure 5. Impeller Hub Blade Loading. 
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Figure 6. Impeller Mean Blade Loading. 
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Figure 7 . Impeller Tip Blade Loading. 
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Figure 8. Vane Oiffuser Entrance Region. 
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Figure 9. Diffuser Cross-section. 
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Figure 10. Leading Edge Region of Vane Diffuser. 
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Figure 11. 90° Annular Bend. 
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Figure 12. Area Distribution for 90° Annular Bend. 
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TABLE IV. 

ANNULAR 

BENO CONTOUR 

COORDINATES. 



R 

Z 

R 

Z 

AREA 

%M/KO 

HUB 

HUB 

SHROUD 

SHROUD 

/ AREA1 

0 

1 4.3000 

5.2 00 5 

14.30 00 

4. 5416 

1.0 000 

10 

1 4.4 538 

5.2130 

14.5600 

4.5450 

1.3404 

20 

14.6137 

5.2 502 

14.3400 

4.5500 

1.1501 

30 

14.7608 

5.3111 

15.03 40 

4. 5720 

1.2 776 

40 

1 4.8966 

5.3943 

15.4190 

4.6580 

1.4524 

50 

1 5.0177 

5.4978 

15.66 70 

4.8520 

1.4911 

60 

15.1211 

5.6189 

15.80 00 

5.1150 

1.3872 

70 

1 5.2C44 

5.7547 

15.8680 

5.4160 

1.2285 

80 

1 5.2653 

5.9018 

15.90 20 

5.69 30 

1.1082 

90 

1 5.3025 

6.0567 

15.9150 

5.9580 

1.0277 

100 

15.3150 

6.2155 

15. 91 84 

6. 2155 

1.0001 


An estimated performance map was prepared for the scaled compressor stage and 
is given in Figure 13. Flow-speed and efficiency lapse rates were maintained 
similar to BU 257 data. However, allowance was made to overall efficiency 
levels to account for the geometrical changes, i.e. reduced area ratio dif- 
fuser and the 90° annular bend. In addition, Reynold's number effects were 
estimated from internal DDA procedures. 

Impeller to shroud clearance distributions for both "build" and "hot running" 
conditions are presented for the 404- III impeller in Figure 14. Build clear- 
ances were deduced from design contours for the "cold" impeller and shroud in 
conjunction with "build" wax check measurements. "Hot running" clearances 
were developed from predicted design speed contours for the impeller and 
shroud and post test rub pin measurements. Scaled values of these clearances 
were assumed in the estimated map shown in Figure 13. 
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Figure 13. Estimated Performance Map for Scaled Compressor. 
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III. STRUCTURAL ANALYSIS 


A complete structural analysis was performed for the scalea impeller. The 
overall analysis was conducted for the defined aerodynamic geometry at a de- 
sign speed of 21789 rpm and with standard day inlet conditions of 14.7 psia 
pressure and 518. 7°R inlet temperature. The material properties were assumed 
to be those of Titanium 6AL4V. The overall structural analysis consisted of 
several individual but complementary tasks, namely: 

o Heat transfer 
o Static stress, and 
o Vibrational analyses. 

The heat transfer analysis was required to provide temperature distributions 
for accurate determination of thermally induced stresses and deflections. 
Boundary conditions consisting of anticipated rig oil temperatures and back 
face seal leakage were obtained from NASA personnel and were included in the 
axisymmetric heat transfer analysis. Results of the heat transfer analysis 
are presented as isotherm lines on the defined scale compressor impeller geom- 
etry in Figure 15. 

The static stress analysis consisted of several tasks: 
o Axisymmetric modeling of the complete wheel geometry 

o Evaluation of stresses from axisymmetric model in terms of low cycle fa- 
tigue and burst margins 

o Triangular plate modeling of the full blade, splitter and backplate 
o Evaluation of triangular plate results in terms of peak stress levels, 
high cycle fatigue and detailed deflection characteristics. 
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Figure 15. Temperature Distributions for Scaled Impeller. 
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Initially, the entire wheel was described with an axisymmetric finite element 
model. This model incorporated the isotherms from the heat transfer analysis, 
the desired wheel geometry (back face contouring and bore diameters) and flat 
radial plates to simulate the weight and approximate stiffness of the impeller, 
blading. 

The axisymmetric finite element model for the scaled impeller geometry is 
shown in Figure 16. Resulting distributions of design speed equivalent, 
axial, radial and tangential stresses are presented in Figures 17, 18, 19 and 
20. The analysis did not_ include an axial clamp load which would result from 
tie bolt stretching. In order to evalaute the possible effects of a given 
axial load on stress levels and distribution, an arbitrarily assumed 10,000 lb 
load (typical of 404-1 II loads) was then applied at the curvic coupling loca- 
tion. The effects of this load are minimal as shown in Figure 21. 

The locations of peak stress were identified and used in a low cycle fatigue 
and burst speed analysis. The results of these analyses are summarized in 
Figure 22 with the highest stresses occurring in the bore. Material proper- 
ties used in the low cycle fatigue analysis are given in Figure 23. As evi- 
dent from the peak stresses of Figure 22 and the properties of Figure 23, 
stress levels were such that lives well in excess of 10^ cycles are pro- 
jected. Burst speeds were calculated to be 182% of design speed based on an 
average radial stress and 199% speed based on average tangential stress. 

Using deflections from the axisymmetric wheel analysis as boundary conditions, 
a triangular plate model was constructed for the full blade, splitter and 
backplate. The full blade model is shown in Figure 24. Figures 25, 26, 27 
and 28 present the results for maximum principal and Von Mises equivalent 
stresses for the full blade pressure and suction surfaces. The same data is 
provided for the splitter in Figures 29 through 33. A segment of the 
backplate was analyzed to evaluate blade/backplate interaction effects. The 
backplate model is presented in Figure 34 with maximum principal and Von Mises 
equivalent stresses given in Figures 35 through 38. 
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Figure 16. Ax i symmetric Stress Model for Scaled Impeller. 


25 




Figure 17. Impeller Equivalent Stresses. 
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Figure 18. Impeller Axial Stresses. 
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Figure 19. Impeller Radial Stresses. 
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Figure 20. Impeller Tangential Stresses. 
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Figure 22. Low Cycle Fatigue and Burst Analysis 
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Figure 24. Triangular Plate Model for Scaled Impeller Full Blade. 
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Figure 25. Pressure Surface ~ Maximum Principal Stress. 
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Figure 26. Pressure Surface ~ Equivalent Stress. 
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Figure 27. Suction Surf ace Maxi munn Principal Stress. 
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Figure 29. Triangular Plate Nbael for Scaled Impeller Splitter Blade. 
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Figure 30. Pressure Surface — Maximum Principal Stress. 
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Figure 31. Pressure Surface Equivalent Stress. 
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Figure 32 . Suction Surface ~ Maximum Principal Stress. 
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Figure 34. 


Triangular Plate Model for Impeller Backplate. 
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Figure 35. Maximum Principal Stress ~ Backf ace Surface. 
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Figure 36. Equivalent Stress ~ Backface Surface. 
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Figure 37. Maximum Principal Stress Flowpath Hub Surface. 
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Figure 38- 


Equivalent Stress ^ 


Flowpath Hub Surface. 
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Maximum equivalent stresses from the triangular plate analysis are summarized 
in Figure 39. Locations A and B are on the full blade leading edge, C is on 
the splitter leading edge and D is at the maximum stress location on both the 
full blade and splitters. Based on the stresses and stress concentration fac- 
tors defined in Figure 39 and the Goodman diagram presented in Figure 40, a 
high cycle fatigue analysis was performed. The maximum steady stress location 
on the splitter yielded the lowest vibratory allowable stress, 14.4 ksi. 

The final output of the static stress analysis was a complete definition of 
the deflection characteristics of the wheel, backplate and blade geometry. 
This definition was used to convert the "hot running" aerodynamic geometry 
into an equivalent manufacturing, "cold", definition. Details of this conver- 
sion and the resulting manufacturing definition are presented in Section IV. 

A summary of pertinent deflections is given in Figure 41. Of special interest 
are points 5 and 2. The calculated axial deflection at the impeller exit was 
0.01463 inches with ground located near the impeller leading edge. This axial 
deflection must be accounted for in the establishment of cold build clear- 
ance. Point 2 is at the location of the curvic coupling. With point 1 
grounded, the change in axial wheel length between points 1 and 2 is 0.00337 
inches with the wheel length shortening due to Poisson's effect. This short- 
ening must be accounted for by initial tie bolt stretch. 

Finite element models similar to those used for static stress analysis were 
prepared for vibrational analysis. Resulting frequencies and mode shapes for 
the first 8 natural frequencies for the full blade and first 3 for the split- 
ter are presented in Figures 42 through 52. A summary of these predicted fre- 
quencies is shown on a frequency-speed diagram for the full blade in Figure 53 
and for the splitter in Figure 54. The first mode of the full blade is well 
above 4th harmonic of rotor speed and is, therefore, not susceptable to inlet 
distortion induced exitation that is often characterized by strong 2nd order 
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Figure 39. Summary of High Cycle Fatigue Analysis. 
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Figure 40. Goodman Diagram for Ti 6AL-4V Material. 
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Figure 41. Scaled Impeller Deflection Summary. 


51 





Figure 42. Full Blade ~ 1st Natural Frequency, 
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Figure 43. Full Blade~2nd Natural Frequency. 
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Figure 44. Full Blaae~3rd Natural Frequency. 


54 













iGLOBfiLiiX^lN) 


Figure 45. Full Blade ^4th Natural Frequency. 
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Figure 46. Full Blaae~'5th Natural Frequency, 
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Figure 47. Full Blade‘~'6th Natural Frequency. 
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Figure 48. Full Blade ~7th Natural Frequency. 
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Figure 49. Full Blade— -8th Natural Frequency. 


59 
















Figure 50. Splitter Blade — ^lst Natural Frequency. 
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Figure 51. Splitter Blaae~2nd Natural Frequency. 
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Figure 52. Splitter Blade ~ 3rd Natural Frequency, 
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Figure 53. Frequency-Speed diagram for Full Blade of Scaled Impeller. 
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Figure 54. Frequency-Speed Diagram for Splitter Blade of Scaled Impeller. 
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(2E0) and its second harmonic, 4E0. The second and eighth modes shown in 
Figure 53 are backplate modes. The remaining modes are essentially blade 
modes with little or no coupling to the backplate structure. 

The sixth diametral mode pattern (60) was selected for analysis of the 
blade/backplate system because it is the only candidate mode which can be 
excited by the 24 engine order (diffuser vane number) when there are 30 air- 
foils on the rotating disk. The first backplate mode (6D-2), Figure 53, is 
acceptable because it is in resonance with 24E0 at a very low speed and should 
be exposed to excitation only during transient operation. The second 
backplate mode (6D-8) should not be of concern unless steady state operation 
near 109% mechanical speed is anticipated. 

In general, the natural frequencies are similar to those of the 404-1 1 1 and, 
therefore, should be totally acceptable for anticipated rig operation. 
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IV. MANUFACTURING DEFINITION OF SCALED IMPELLER 


This section describes the procedure by which the "hot running" impeller is 
converted into the "cold" or manufacturing definition. The initial portion of 
this process is to combine the "hot running" impeller geometry defined by 
aerodynamics with the deflection characters itics calculated in the static 
stress analysis. 

The "hot to cold" static stress analysis is performed by iteratively adjusting 
the impeller geometry ano analytically "spinning" this geometry to design 
speed and temperature. Convergence is assumed when this "hot running" geom- 
etry matches that originally defined by the aerodynamics group. The final 
output of this analysis is a full 3-D cartesian definition which gives the X, 
Y, Z location of the desired "hot running" blade and the V X, V Y, 7Z deflec- 
tion required to convert this "hot running" blade to the manufacturing defini- 
tion. 

Blade generation programs, then, apply these deflections, reconstruct the 
"cold" blade, modify the shroud contour for the desired "hot running" clear- 
ance and, finally, define the cold blade on a series of planes passed both 
parallel ano perpendicular to the engine axis. 

Figure 55 indicates DDA practice for applying the "hot running" clearance. A 
constant clearance is first subtracted from the "hot" shroud line defined pre- 
viously in Table II. The original 404- III had a constant 0.010 inches removed 
from the impel ler^s^roud line. Using the 1 .6529 linear scale factor, a clear- 
ance of .0l6529*was used for the scaled impeller. Since a static structure 
analysis was not performed, an absolute definition of the required "cold" 
build clearance is not possible. However, the change in clearance (relative 
to the originally defined "hot" shroud line) calculated for the impeller is 
provided in Figure 56 for reference. 

^ To .00%" Her CLeM 7/IX/23 
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Figure 55. OOA Procedure for Clearance Allowance. 
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Figure 56. Scaled Impeller Clearance Change from "cold build" to 


"hot running" Condition. 
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The final "cola, manufacturing" definition of the impeller blaae and splitter 
with the 0.016529 inch clearance removed is given in Tables V and VI, respec- 
tively. The "cold" blade surfaces are assumed to be constructed by straight 
lines from the hub to the shroud along the defined quasi-normals. The split- 
ter blades are equally spaced between full blades at the impeller exit. 

Planar sections were then passed through this "cold" geometry. The printed 
coordinate intersections of these planes with the quasi-normals is called a 
blade book and is included as Attachment A. Punched card definition of the 
hot and cold blade surfaces defined in Tables II, III, V and VI and the blade 
book coordinates are included as Attachment B. Mi 1 ar plots of the planar 
intersections are know as Master Charts and are included as Attachment C. 
These plots are five times scale and can be used for final part inspection. 
Finally, a SK drawing was prepared to define the wheel geometry ana locations 
of the planar Master Chart sections. This drawing is Attachment 0. 
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